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ABSTRACT 

Both the Fe II UV emission in the 2000-3000 A region (Fe II(UV)) and 
resonance emission line complex of Mg II at 2800 A are prominent features in 
quasar spectra. The observed Fe II(UV)/Mg II emission ratios have been pro- 
posed as means to measure the buildup of the Fe abundance relative to that 
of the a-elements C, N, O, Ne and Mg as a function of redshift. The current 
observed ratios show large scatter and no obvious dependence on redshift. Thus, 
it remains unresolved whether a dependence on redshift exists and whether the 
observed Fe II(UV) /Mg II ratios represent a real nucleosynthesis diagnostic. We 
have used our new 830-level model atom for Fe + in photoionization calcula- 
tions, reproducing the physical conditions in the broad line regions of quasars. 
This modeling reveals that interpretations of high values of Fe II(UV)/Mg II are 
sensitive not only to Fe and Mg abundance, but also to other factors such as 
microturbulence, density, and properties of the radiation field. We find that the 
Fe II(UV)/Mg II ratio combined with Fe II (UV)/Fe II (Optical) emission ratio, 
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where Fe II (Optical) denotes Fe II emission in 4000-6000 A band, can be used 
as a reliable nucleosynthesis diagnostic for the Fe/Mg abundance ratios for the 
physical conditions relevant to the broad-line regions of quasars. This has ex- 
treme importance for quasar observations with the Hubble Space Telescope and 
also with the future James Webb Space Telescope. 

Subject headings: atomic processes — line: formation — methods: numerical- 
quasars: emission lines 

1. INTRODUCTION 

The observed ratio of the ultraviolet Fe II UV emission flux (2000—3000 A; hereafter 
Fe II (UV)) to that of the Mg II resonance doublet at 2800 A in broad- line regions (BLRs) 
of quasars has a strong potential for being a fundamental cosmological metallicity indicator. 
In models of galactic chemical evolution, the ratio of Fe to a-element (O, Ne, Mg) abun- 
dances constrains the ages of star-forming systems (Wheeler et al. 1989). The Fe/ 'a age 
constraint follows from a different enrichment timescale, namely: a-elements are produced 
in supernova explosions of short-lived massive stars (primarily SN Type II), while Fe has 
a large contribution from the longer-lived intermediate-mass binaries that produce Type la 
supernovae. Predictions indicate that the Fe/a abundance ratio should increase by more 
than an order of magnitude for stellar systems older than 1 Gyr (Hamann & Ferland 1993; 
Yoshii et al. 1996). However, if star formation has a significant contribution from very 
massive objects (100 < M* < 1000M Q ) as suggested recently (Heger & Woosley 2002), then 
substantial Fe production can occur at times earlier than 1 Gyr. The answer to whether the 
Fe/a abundance ratio has a sharp break at 1 Gyr has important consequences not only for 
the evolution of the elements, but possibly for understanding the formation of active galactic 
nuclei (AGN) as well. Whether the Fe/a enrichment is due to SN la or SN II explosions is 
a key problem to be addressed in constructing our picture of galactic evolution in the Early 
Universe. 

The delay time in Fe enrichment can be estimated from observations by plotting the 
observed Fe II (UV) /Mg II emission ratio (or the Fe/Mg abundance) as a function of redshift. 
Recent Fe II (UV)/Mg II emission ratio measurements (e.g. Iwamuro et al. 2002; Dietrich et 
al. 2002; Freudling et al. 2003) show a wide range of values. That obscures any dependence 
on redshift. It is still unclear: (1) whether a dependence on redshift exists; and (2) whether 
the derived observed Fe II/Mg II values are sensitive to other physical parameters besides 
Fe abundance. 
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How the observed Fe II (UV) /Mg II emission ratio varies with physical conditions must 
be resolved before one can develop a reliable model that accurately predicts the Fe II emission 
line spectra in the BLRs of luminous AGNs. The large velocity dispersion of Fe II emission in 
BLRs, the large number of overlapping lines, and the overall richness of the Fe II spectrum 
(Verner et al. 1999) form broad emission blends, producing a "Fe II pseudo-continuum" 
and/or prominent features from the UV through visual (4000—6000 A; hereafter Fe II (Op- 
tical)) to IR. This necessitates a complete simulation of the physical processes affecting the 
Fe + ion to deduce how the Fe II emission varies with properties of the radiation field, density, 
temperature, and Fe abundance in the emitting regions. To investigate what really makes 
Fe II emission strong in the rest-UV, we have used a new Fe + model ion. The model includes 
830 atomic levels for Fe + (up to ~ 14 eV) and represents an expansion of the previous model 
of Verner et al. (1999). 

In this paper we evaluate how Fe II(UV)/Mg II and Fe II(UV)/Fe Il(Optical) emission 
ratios vary due to extreme changes of the radiation field. Then, for the given density and 
radiation field, we investigate the Fe II emission properties within a reasonable range of Fe 
abundance and microturbulence. 



2. EFFECTS OF PHYSICAL CONDITIONS ON Fe II EMISSION 

2.1. The expanded Fe + ion model 

Our earlier 371-level Fe + ion model (for energies below 11.6 eV), incorporated into the 
code CLOUDY 8 , predicts the emission fluxes of 68,635 Fe II lines (Verner et al. 1999). 
This 371-level model, plus accurate knowledge of the radiation field in a low-density plasma 
successfully predicted Fe II emission in the lower density environments of the Orion Nebula 
and Eta Carinae Ejecta (Verner et al. 2000; Verner et al. 2002). Comparison of observations 
with model results indicates that anomalously strong Fe II emission lines can be the result 
from line pumping due to coincidence of Fe II transitions with strong emission lines of other 
abundant ions (H I, He II, Si III, C IV, etc.; Hartman & Johansson 2000). 

The higher radiation and particle density regimes in the BLRs of quasars favor the pop- 
ulation of higher energy levels in Fe + . Consequently, we expanded the original Fe + model 
to 830 levels (energies below 14.06 eV yeilding 344,035 transitions) and included more accu- 
rate atomic data - new experimental energy levels (Johansson 2002, private communication) 
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and A- values from the Amsterdam database 9 (see also Raassen & Uylings 1998). Where 
A-values were absent for many semi-forbidden transitions, and rate coefficients uncertain for 
Lya pumping channels and primary cascades, we have inserted reasonable estimates. 

The complete simulation of the physical processes affecting the Fe + ion makes it possible 
to deduce the density, temperature, plus predict iron abundance in the emitting regions. To 
accurately model the Fe II emission and Fe abundances, we must quantify the effects of the 
radiation field, and, only then, use Fe II(UV)/Mg II ratios to determine relative abundances. 
This approach is especially important for high redshift quasars where selection effects due to 
the luminosity may play a significant role. Thus, our knowledge about radiative excitation 
in Fe + is crucial to predict the Fe II spectrum and Fe abundances at high redshifts. The 
questions is, "Can we use the observed Fe II(UV)/Mg II and Fe II(UV)/Fe Il(Optical) ratios 
to derive Fe abundances?" 



2.2. Model parameters and comparison of the 830 and 371-level models 

In our modeling, we have chosen physical conditions where lines from low-ionization 
species like Fe II dominate the emission. We adopt a range of Fe abundance and a constant 
density, n H = 10 9,5 cm -3 , with a total hydrogen column density, N H = 10 24 cm -2 . We further 
assume that the flux of hydrogen ionizing photons at the illuminated face is 10 20 5 cm -2 s -1 . 
These parameters for BLR conditions are within the range of values taken from Verner et 
al. (1999), and Verner (2000). We employ the characteristic AGN continuum described in 
Korista et al. (1997), which consists of a UV bump peaking near 44 eV, a f u oc u^ 1 X-ray 
power law, and a UV to optical spectral index, a ox = —1.4. For the adopted parameters, 
we have investigated how variations in abundance and microturbulence alter the intrinsic 
emission ratios of Fe II(UV)/Mg II and Fe II (UV)/Fe Il(Optical). 

The differences in the predicted Fe II spectrum calculated using 371 levels and 830 
levels are negligible for lower densities and radiation field with depleted iron abundances 
(e.g. Orion Nebula). Differences in the predicted spectra from the two models become 
more significant with increased radiation field and iron abundance. Moreover, the shape of 
the radiation field in AGNs leads to increased population of higher energy levels in Fe II. 
Since densities are higher and the radiation field is more intense in BLRs, these upper levels 
provide additional excitation channels in the far UV (energies higher than 12 eV). 

Figure 1 compares the predictions using our 371- and 830-level models for the Fe II 
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pseudo-continuum appropriate for BLR physical conditions in QSOs. The plotted spectra are 
given as flux versus wavelength, where flux is uf u , in units of erg cm~ 2 s _1 . The comparison 
clearly shows how the increased number of atomic levels produces a stronger Fe II pseudo- 
continuum in QSO spectra. The 830-level ion model is far more accurate than the previous 
best efforts (Wills et al. 1985; Verner et al. 1999). The increased number of Fe + energy 
levels influences not only the Fe II spectrum but lines of other elements and the whole energy 
budget also. However, the complete study of the new theoretical emission spectrum is not a 
subject of this paper. 

Table 1 shows the predicted Fe Il(UV)/Mg II and Fe II(UV)/ Fe Il(Optical) ratios for 
830-level and 371-level models of Fe II versus contributions from the radiation field for a 
solar iron abundance and microturbulence of 1 km s _1 . The models were calculated for 
radiative continuum pumping and for pure collisional excitation. There is almost no change 
in predictions of 371-level and 830-level models for pure collisional excitation, if pumping 
by the radiation field is ignored. However, the predicted ratios are strongly influenced by 
the radiation field and are sensitive to contributions from the upper levels in Fe II. Figure 2 
demonstrates the predicted Fe II spectrum due to collisions and radiative pumping. 



2.3. Effects of microturbulence and abundance 

How does the Fe II emission depend upon microturbulence? The first successful models 
were developed by Netzer & Wills (1983) who also demonstrated that the Fe II emission 
strength is sensitive to the turbulent velocitiy. This sensitivity was confirmed by modeling 
by Verner et al. (1999), who showed that the predicted UV Fe II emission indeed became 
stronger with increased microturbulence. To explore how microturbulence affects the ob- 
served Fe II(UV)/Mg II measurements, we calculated models with different microturbulence 
parameters (v tur b = 1, 5, 10, 100 km s _1 ). 

Figure 3 shows Fe II(UV)/Mg II and Fe II(UV)/Fe Il(Optical) ratios as a function of mi- 
croturbulence. For small values of v tur b and with solar abundance, the Fe II(UV) /Mg II ratio 
is smaller than that of Fe II(UV) /Fe Il(Optical). The ratios are the same at Vturb = 5 km s _1 . 
However, the Fe Il(UV)/Mg II ratio grows faster than that of Fe II(UV)/Fe Il(Optical) as a 
function of microturbulence. The UV Fe II transitions often arise from levels at high exci- 
tation energies, and the density of levels increases substantially above 10 eV. Thus, the UV 
transitions are closer in energy than those in the optical range, and their strengths are more 
affected by microturbulence than those of optical lines. Although both the Fe II(UV)/Mg II 
and Fe II(UV)/Fe Il(Optical) ratios increase with microturbulence, dependencies are not 
identical. Our calculations demonstrate that the most reasonable values for Fe II(UV)/Mg II 
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ratios correspond to v tur b = 5 — 10 km s l . 

Is there any specific Fe II feature, line or spectral index in AGN spectra that can be 
used to derive the iron abundance accurately? To probe the effect of varying abundance, we 
repeated the calculation with one and five times solar iron abundance, keeping the micro- 
turbulence constant at v tur b = 1 km s -1 . We define the iron abundance sensitivity ratio 
as (I(Fe II) 5 - I(Fe Il)i)/I(Fe II)! and signify integrated fluxes I(Fe II)! and I(Fe II) 5 for 
abundances of solar and five times solar. Calculations with v tur b — 10 km s _1 (Figure 4) 
show that this sensitivity ratio has the same slope as it has with v tU rb = 1 km s -1 . 

We find that the UV Fe II lines are generally less sensitive to abundance than those 
in the optical range. Energy levels responsible for UV and optical lines are not populated 
proportionally with increasing abundance. The structure of the Fe + ion favors the production 
of more optical emission through cascades from upper levels. 

Figure 5 shows how Fe II(UV)/Mg II and Fe II(UV)/Fe Il(Optical) emission ratios 
change with iron abundance (1, 5, and 10 solar). The increase in Fe II (Optical) emis- 
sion with the Fe abundance is stronger than that in Fe II (UV) emission. Comparisons 
of radiation and/or turbulence effects confirm that the abundance is not the only factor 
that affects the Fe II (UV)/Mg II ratio. A factor of 5 in abundance increases the Fe II 
(UV)/Mg II ratio by factor less than 2. By using Fe II (UV)/Mg II in combination with 
the Fe II(UV)/Fe Il(Optical) emission ratio, we can deduce abundances while accounting 
for microturbulence and radiation field effects. Table 2 shows predictions for models with 
different microturbulence (v tur b = 1 — 100 km s _1 ) and iron abundance (1 — 10 of solar). The 
columns are the following: microturbulence, abundance, and calculated Fe II(UV)/Mg II and 
Fe II(UV)/ Fe Il(Optical) ratios. 



3. SUMMARY 

By applying the new 830- level model ion for Fe + , which is incorporated into detailed 
photoionization calculations, we have probed the feasibility of using the strong Fe II emission 
seen in quasar spectra to derive accurate Fe/Mg abundance enrichments. This Fe II emission 
would be especially important in determining the enrichment as a function of redshift and 
the nature of star formation in the Early Universe. We conclude that abundance is not the 
only factor that makes Fe II emission strong over a wide range of wavelengths and is not even 
the strongest factor. Our modeling indicates that microturbulence also leads to preferential 
strengthening of the Fe II emission in UV. 

On the basis of the calculations presented here, we make five basic conclusions: 
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1. In contrast to recent observational Fe II/Mg II measurements, we emphasize the need 
for low-z quasar observations to calibrate detailed modeling as used here. Our model clearly 
demonstrates that a simple Fe II template cannot be applied for all quasars, because the Fe 
II sensitivity to radiation and turbulence is much stronger than that of the Mg II doublet. 
Only detailed modeling that includes the all these effects can provide more accurate Fe II 
templates. 

2. Due to strong Fe II line overlap, we must adopt a wavelength band approach to 
get constraints on density, radiation field, and turbulence in BLRs. This approach is more 
complex, but the computational tools are readily available. 

3. Further modeling, beyond that presented here, is necessary. Model comparisons with 
high S/N observations of quasars, and of even higher resolution observations of objects like 77 
Car, will produce better atomic parameters, and more detailed understanding of the effects 
of the radiation field in astronomical objects. Future work using improved atomic data to 
fully explore effects of the radiation field is in progress. 

4. These results clearly illustrate that the continuum of AGNs is greatly affected by 
Fe II. The spectrum of Fe II simply cannot be ignored. Understanding Fe II is crucial for 
obtaining realistic modeling of restframe AGN spectra from the UV through the IR. Studies 
in the expanded spectral range, much wider than used in any previous efforts, are required. 

5. We suggest that our model is to be applied to low-redshift quasars spanning both a 
wide luminosity range, and the 1,000 — 10,000 A restframe wavelength range. As the first 
step in this direction, the STIS and NICMOS capabilities can be used immediately to study 
low-redshift quasars. The Near Infrared Spectrograph (NIRSpec) on board the James Webb 
Space Telescope (JWST) would make it possible to measure spectra from 0.6 to 5 micron 
with resolving power in the range needed to obtain Fe II(UV) and Fe Il(Optical) bands from 
low-redshift to high-redshift quasars. 
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Fig. 1. — Fe II pseudo-continuum predicted by 371 (inset) and 830 level models for Fe + in 
QSOs' BLR conditions (section 3.1). The 1000—7000 A wavelength range devided into small 
bins of 6 A (consistent with a 500 km s™ 1 FWHM at 3600 A) and Fe II flux calculated 
in each bin. The total Fe II flux in range 2000—3000 A increased by factor 1.4 (Table 1) 
compared to the older model (Verner et al. 1999). The dashed line separates the Fe II 
pseudo-continuum from strong lines in the 371-level model. The 830-level model predicts 
many strong Fe II features in the 1600—3500 A range. 

Fig. 2. — Fe II spectrum calculated in a radiative continuum pumping model (solid line) and 
a pure collisional model (dotted line and inset). Broadening by a non-turbulent FWHM of 
2,000 km s^ 1 (as observed for broad lines in the quasar 3C 273) has been used to enable 
comparison with observations. 

Fig. 3.— Fe Il(UV)/Mg II and Fe II(UV)/Fe Il(Optical) ratios as a function of microturbu- 
lence. See text for further discussion. 

Fig. 4. — The general sensitivity of the Fe II continuum to changes in abundance (for 
constant v tur b = 10 km s -1 ). The Fe II(UV) emission is less sensitive to abundance than 
Fe Il(Optical). 

Fig. 5. — Fe II(UV)/Mg II and Fe II(UV)/Fe Il(Optical) ratios as a function of abundance. 
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Table 1. Predicted Fe II(UV)/Mg II & Fe II(UV)/Fe II (Optical) Emission Ratios versus 
Excitation for Solar Abundance and Turbulence of 1 km s" 1 . 



Number of levels Excitation Fe II(UV)/Mg II Fe II(UV)/Fe Il(Optical) 

371 collisional 0.5 2.8 

radiative 4 11.4 

830 collisional 0.3 2.8 

radiative 5.4 8 
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Table 2. Predicted Fe II(UV)/Mg II & Fe II(UV)/Fe II (Optical) Emission Ratios versus 

Abundance and Turbulence. 

Turbulence (km s- 1 ) Abundance Fe II(UV)/Mg II Fe II(UV)/Fe II (Optical) 



1 1 5.4 8 

1 5 9.3 6.5 

1 10 10.7 5.6 

5 1 10 10 

10 1 13.5 10.7 

100 1 30 18 

5 5 14 6.8 

10 5 20 7.3 

10 10 27.5 6.4 

100 10 77 9.2 



